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ABSTRACT 

We report the discovery of pulsed 7-ray emission and X-ray emission from the 
black widow millisecond pulsar PSR J2051-0827 by using the data from the Large 
Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope and the 
Advanced CCD Imaging Spectrometer array (ACIS-S) on the Chandra X-ray 
Observatory. Using 3 years of LAT data, PSR J2051-0827 is clearly detected 
in 7-ray with a signicance of ~ 8cr in the 0.2 - 20 GeV band. The 200 MeV - 
20 GeV 7-ray spectrum of PSR J2051-0827 can be modeled by a simple power- 
law with a photon index of 2.46 ± 0.15. Significant (~ 5<x) 7-ray pulsations 
at the radio period were detected. PSR J205 1-0827 was also detected in soft 
(0.3-7 keV) X-ray with Chandra. By comparing the observed 7-rays and X-rays 
with theoretical models, we suggest that the 7-ray emission is from the outer 
gap while the X-rays can be from intra-binary shock and pulsar magnetospheric 
synchrotron emissions. 

Subject headings: gamma rays: stars — (stars:) pulsars: individual (PSR J2051- 
0827)— X-rays: individual (PSR J2051-0827) 



1. INTRODUCTION 

It is widely accepted that millisecond pulsars (MSPs) have been spun up by accreting 
masses from a close binary companion (Alpar et al. 1982). A sub-class of binary systems 
consist of a MSP and a very low mass companion (M c < 0.05M SU „) orbiting in a tight orbit 
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(P b < 24/ir). In such systems particles or 7-rays emitted from the pulsar are thought to 
be ablating the companion, which would eventually destroy the companion, resulting in a 
solitary MSP (Fruchter et al. 1995; Takata et al. 2011); MSPs in this kind of binary system 
are usually called Black Widow pulsars. 

PSR J2051-0827 is the second black widow pulsar after PSR B 1957+20 (Fruchter, Stine- 
bring & Taylor 1988). PSR J2051-0827 was discovered in 1996 in the Parkes all-sky survey of 
the southern sky for low-luminosity and millisecond pulsars (Stappers et al. 1996). The spin 
period of PSR J2051-0827 is 4.5ms. The MSP is orbiting around a very low mass companion 
(m c < 0.027M SU „) with a orbital period of 2.4 hrs, thus joining the black widow populations 
(Stappers et al. 1996). Since the discovery of this pulsar, follow-up observations have been 
performed to reveal the nature of this black widow system (Stappers, Bessell & Bailes 1996; 
Stappers wt al. 1999; Stappers et al. 2001b; Doroshenko et al. 2001; Lazaridis et al. 2011). 

Before the launch of the Fermi-LAT , all black widow pulsars have been found by radio 
surveys but most of them are in globular clusters (GCs). As GCs usually hold a large 
number of MSPs and many GCs are 7-ray emitters (Abdo et al. 2010; Kong et al. 2010; 
Tam et al. 2011). However, we cannot resolve individual pulsars in 7-rays in GCs due 
to the limited angular resolution of Fermi-LAT , which makes studies of 7-ray black widow 
pulsars in GCs challenging. After the launch of the Fermi-LAT, more and more field black 
widows were identified as Fermi sources, opening a new observation window of Black Widow 
pulsars studies, providing a larger sample for characteristic studies on black widow pulsars. 
So far a total of 10 black widow pulsars have been found in the Galactic field, in which 7 of 
them are Fermi sources and at least 3 of them show pulsed 7-rays (Roberts 2011). Among 
those which has remained undetectable in 7-rays, PSR J205 1-0827 is the closest one with a 
dNE2ooi distances of 1.04 kpc (Lazaridis et al. 2011), which increases the chance of detecting 
it in 7-rays despise it has a relativity low spin down power E = 0.53 x 10 34 erg/s (Stappers 
et al. 1996). In this letter, we present the discovery in 7-ray and X-ray using Fermi-LAT 
and Chandra data. 



2. DATA ANALYSIS AND RESULTS 

2.1. GAMMA-RAY DATA 

In our analysis, we used the LAT data between 2008 August and 2011 September (~3 
years of data). To reduce and analyze the data, the Fermi Science Tools v9r23pl package, 



Fig. 1— 2°x2° LAT count map in 1-20 GeV centered at PSR J2051-0827 with a pixel 
size of 0.025°, smoothed by a Gaussian kernel of 0.15°. The cross in magenta represents 
the nominal position of PSR J2051-0827. A scale bar is used to indicate counts/pixel. No 
background subtraction has been applied. 

which is available from the Fermi Science Support CenteiQ was used. We used Pass 7 data 
and selected events in the "Diffuse" class (i.e. event class 2) only. In addition, we excluded 
the events with zenith angles larger than 100° to reduce the contamination by Earth albedo 
gamma-rays. The instrumental response functions (IRFs) "P7SOURCE_V6" were adopted 
throughout this study. 

Events were selected within a circular region-of-interest (ROI) centered at the nominal 
position of PSR J2051-0827 (i.e. RA=20 /l 51 m 07 s .51 Dec=-08°27'37".7). In order to reduce 
systematic uncertainties and achieve a better background modeling, a circular ROI with a 
diameter of 10° was adopted throughout the unbinned likelihood analysis. For a preliminary 
inspection of the chosen part of the sky around PSR J2051-0827, a binned photon count 
map in 1 — 20 GeV was produced with task gtbin (see Figure [T]). 7— ray excess can be 
clearly identified around the nominal position of PSR J2051-0827 even before subtracting 
the background. 

To investigate the spectral characteristics of PSR J2051-0827, we performed an un- 
burnned likelihood analysis with gtlike, by assuming a Power Law (PL) as well as Power 
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Law with an exponential cutoff (PLE) model for a point 7-ray source at the nominal po- 
sition of PSR J2051-0827. Photon energies are restricted to 200 MeV to 20 GeV. In the 
background model, we included the Galactic diffuse model (gal_2yearp7v6_v0 . f its) and 
the isotropic background (iso_p7v6source.txt), as well as all point sources reported in 
the Second Fermi-LAT Source Catalog (2FGL) within 10° from the center of the ROI. All 
these 2FGL sources were assumed to be point sources which have the specific spectrum 
suggested in the 2FGL catalog (Abdo et al. 2011). While the spectral parameters of the 
2FGL sources locate within the ROI were set to be free, we kept the parameters for those 
lie outside our adopted ROI fixed at the values given in the 2FGL catalog (Abdo et al. 
2011). We also allowed the normalization factors of the two diffuse background components 
to vary. The best-fit PL model yields a photon index of T — 2.46 ± 0.15 and a test-statistic 
(TS) value of 66 which corresponds to a significance of ~8<r. Its photon flux in this band 
is found to be (6.62 ± 1.29) x 10 -9 cm -2 s _1 The corresponding integrated energy flux is 
/ 7 = 5. 92^3^4 x 10~ 12 erg cm' 2 s _1 F For the PLE model, the best fit results in a photon 
index of T = 1.75 ± 0.58 and a cutoff energy of E cuto ff = 2.5 ± 2 GeV, giving a TS value 
of 70 which corresponds to a significance of ~8.4<r. Based on the likelihood ratio test, the 
PLE model are not statistically required to fit the observed 7-ray spectrum. This may be 
due to the low photon statistic and the contamination from the diffuse 7-ray background. 
The binned spectrum is displayed in Figure [2} 

As PSR J205 1-0827 is not reported in the 2FGL catalog (Adbo et al. 2011), we tried 
to reproduce the non-detection result using binned likelihood analysis by defining a 7° x 7° 
square ROI, together with the same time span (2 years) and energy cut (0.1-300 GeV) as in 
the 2FGL catalog (Abdo et al. 2011). The resulting TS is below 25, lower than the selection 
criteria set by the 2FGL catalog (TS > 25). Therefore it is expected that PSR J2051-0827 
was not in the 2FGL catalog. After reproducing the non-detection result, we expand the 
time span by using ~3 years of photon data to perform a binned likelihood analysis again 
to check whether PSR J2051-0827 will be detectable with a longer exposure. By assuming 
a simple power law model, the best-fit model yields a photon index of T — 2.39 ± 0.13 
and a test-statistic (TS) value of 49. The photon index is consistent with the unbinned 
likelihood analysis using 2 years of photon data with the same energy cut, and therefore 
PSR J2051-0827 is detected with ~3-years of LAT data. 

We generated a 2°x 2° TS map in 1-20 GeV centered at the nominal position of PSR 
J2051-0827 by using gttsmap. This is shown in Figure [3] With the aid of gtfindsrc, We deter- 
mined the best-fit position of PSR J2051-0827 in 1 - 20 GeV to be RA=20 /l 51 m 17.28 s Dec=- 



°Thc quoted errors of the energy flux have taken the statistical uncertainties of both photon index and 
prefactor into account. 



Fig. 2.— Fermi LAT spectrum of PSR J2051-0827 in the 0.2 - 10 GeV band. Curve and 
straight line represent the best-fit PLE & PL model respectively. 

08°26 22.6 (J2000) with la error radius (statistical) of 0.09° which is illustrated with a black 
circle in Figure [3| 

By using the latest radio timing ephemeris reported from Lazaridis et al. (2011), we 
searched for 7-ray pulsations on PSR J2051-0827 by phase binning Fermi LAT photons 
received within ~3 years photon data using the Fermi plugin for TEMP02, Fermi LAT 
photons were corrected to the solar system barycenter using the JPL DE405 solar system 
ephemeris. To achieve the best 7-ray pulse profile, we test on two parameters: the radius 
of the aperture (0.1° < r < 1°) and the lowest photon energy (100 < E min < 1000 MeV), 
with an E max at 300 GeV. The best signal for 7-ray pulsation of PSR J2051-0827 (r = 
0.5°, E min > 200 MeV) is presented in Figure [1] with a H-test TS of 35 (de Jager & Biisching 
2010), corresponding to a detection significance of 4.9 a. We are therefore confident that the 
newly discovered 7-ray source is associated with PSR J2051-0827. 

Since orbital modulations was observed in both radio and optical band, we attempted 
to search for any evidence on orbital modulations by phase-folding the 7-ray photons with 
the orbital parameter using TEMP02 but no positive result has been found using the orbital 
parameter obtained by Lazaridis et al. (2011). 

2.2. X-RAY DATA 

PSR J2051-0827 was observed 5 times with the Chandra X-ray Observatory on 2009 
March 22, 30, and July 5 (ObsID 10106-10110). All the observations were taken using the 



Fig. 3. — Test-statistic (TS) map in 1-20 GeV of a region of 2°x2° centered at the nominal 
position of PSR J2051-0827 (magenta cross). The color scale that used to indicated the TS 
values is shown by the scale bar below. The circle in black represents the la positional error 
circle determined by gtfindsrc. 

Advanced CCD Imaging Spectrometer array (ACIS-S) with the very faint mode. Data were 
collected with a frame transfer time of 3.2 sec and the exposure time for each observation 
is about 9 ksec. We used CIAC0 version 4.3 and XSPECQ version 12.7 packages to perform 
data reduction and analysis. We reprocessed the raw data to apply the most up-to-date 
calibration and to make use the very faint mode. In order to reduce the background, we 
restricted the photon energies between 0.3 and 7.0 keV for all our data analysis. We also 
inspected the background count rates from the SI chip and no flaring event was found in all 
the observations. 

PSR J205 1-0827 is marginally seen in each observation with less than 10 counts. We 
therefore combined all 5 observations with a total exposure time of ~ 44 ksec for subsequent 
analysis. We extracted the energy spectrum from a 2 arcsec circular region centered on 
PSR J2051-0827; a total of 43 counts were extracted. For the background, we selected a 
source-free region with a radius of 30 arcsec, which results 1.1 background count in the 
source region. We therefore did not subtract the background for spectral analysis. Response 



6 http://cxc.harvard.edu/ciao/ 
7 http://heas arc. nasa.gov/xanadu/xspec/ 
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Fig. 4.— 7-ray light curve of PSR J2051-0827 using Fermi plugin for TEMP02 . 



matrices of each observation were generated and co-added with CIAO. Because of the low 
count statistic, we fit the unbinned spectrum with an absorbed power-law and an absorbed 
blackbody models. We fixed the absorption at the Galactic value (Nh = 6 x 10 20 cm -2 ) and 
used CASH statistics (Cash 1979) to estimate the best-fit parameters and their errors. Both 
models provide an acceptable fit. The best-fit photon index is 2.51 ± 0.47 (90% confidence) 
yielding an unabsorbed 0.3-10 keV flux of 8.5 x 10~ 15 ergs cm -2 s -1 . For the blackbody 
model, the best-fit temperature is kT = 0.25 ±0.05 keV and the unabsorbed 0.3-10 keV flux 
is 4.9 x 10~ 15 ergs cm -2 s -1 . We also binned the data with at least 5 counts per spectral 
bin and used x 2 statistics to derive the best-fit model. The results are consistent with those 
using CASH statistics. 



3. DISCUSSION 

We report the discovery of the 7-ray pulsation and the X-ray detection from the second 
black widow system PSR J2051-0827 using Fermi-LAT and Chandra-AClS-S data. We 
compute the rotation energy loss, E = 4n 2 IP/P 3 , where I is the pulsar's moment of inertia. 
We assume / = 10 45 g cm 2 for an ordinary MSP, which gives a E — 5.49 x 10 33 erg s _1 . 
By adopting a distance of 1.04 kpc from Lazaridis et al. (2011), the 7-ray and the X-ray 
luminosity are found to be L 7 = 7.66 x 10 32 erg s~ x and Lx = 1.01 x 10 30 erg s _1 in 0.2-20 
GeV and 0.3-10 keV respectively. Therefore the 7-ray luminosity is ~14% of the spin down 



-8 - 



power of PSR J2051-0827, consistent with other MSPs in the Galactic field. 

Comparing PSR J2051-0827 with those black widow pulsars with both 7-ray and X-ray 
emission detected [PSR B1957+20 (Huang & Becker 2007; Huang et al. in prep.), PSR 
J2214+3000 (Ransom et al. 2011) and PSR J2241-5236 (Keith et al. 2011)], they all show 
similar X-ray properties in terms of their black body temperature (~ 0.25 keV) and X-ray 
luminosity (~ io 30 ~ 31 erg s _1 ). It is worth noting that except for PSR B1957+20, all sources 
have only less than 90 X-ray photons for spectral analysis and blackbody model may not 
be the only acceptable spectral component. For instance, we cannot rule out a power-law 
model for PSR J205 1-0827. Furthermore, a deep Chandra observation of PSR B1957+20 
shows that an additional power-law component is required (Huang et al. in prep.). Unlike 
the other three black widow systems, PSR J2051-0827 is a previously uncatalogued 7-ray 
source and the discovery of the 7-ray emission is entirely driven by the radio observations. 
In all four 7-ray and X-ray emitting black widow pulsars, X-ray pulsation was only detected 
for PSR B1957+20 (Guillemot et al. 2011). 

Takata et al. (2011) study the X-ray and 7-ray emissions from the millisecond pulsar 
in black widow systems. They discuss 7-ray emission from the outer gap accelerator and 
predict its luminosity to be L 7 ~ 8 x 10 32 ergs -1 for PSR J2051-0827, which is consistent 
with the results of the present observation. The observed X-ray emissions from black widow 
system must be composed of several components. For example, Guillemot et al. (2011) 
found that the first 7-ray black widow system, PSR B 1957+20, shows that the total X-ray 
emissions are composed of the non-pulsed (~ 67%) and pulsed emissions (~ 33%). The 
non-pulsed emission will be originated from the synchrotron emissions from the particles 
accelerated at the intra-binary shock due to the interaction between the pulsar wind and the 
matter injected from the companion star. For the emissions from the pulsars, the emissions 
from the heated polar cap region due to the return current or from the magnetospheric 
synchrotron emissions, or both components can contribute to the observed emissions. 

If the X-ray emissions are from an intra-binary shock, Takata et al. (2011) estimated 
L x (0.5 - lOkeV) ~ 2 x 10 30 erg/s for PSR J2051-0827, if the magnetization parameters (the 
ratio of the magnetic energy and the particle energy), is a ~ 0.1 and the power index of 
energy distribution of the emitting particles is ~ 3. Such a model will produce an X-ray 
spectrum with a photon index of ~ 2 which is consistent with the observed value. For the 
magnetospheric synchrotron emissions, Takata et al. (2011) estimate L x ~ 5 x 10 29 erg/s, 
which may be slightly smaller than the observation. For 0.3-10 keV bands, so called core 
component of the heated polar cap emissions, which has a temperature ~ 10 6 K and the 
effective radius 10 4 cm, will contribute to the observed emissions. Using equation (28) in 
Takata et al. (2011), the temperature and luminosity of the core component for PSR J2051- 



-9- 



0827 are estimated as ~ 0.25 keV and L x ~ 4 x 10 30 erg/s, respectively, where the effective 
radius is assumed to be R ~ 0.1 km. This is consistent with the observed values. 

Note that Lazaridis et al. (2011) have shown the secular variations on orbital period on 
PSR J2051-0827, in which the timing model is valid until MJD 54888, which do not cover 
the whole range of the Fermi-LAT data used in folding the 7-ray light curve. Although 
there is no indication of any orbital period variation in the 7-ray light curve but as the 
orbital parameter are subjected to change over time, further radio/optical observations are 
necessary to refine the 7-ray pulse profile. 
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